The thermopower as a fingerprint of the Kondo breakdown quantum critical point 
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We propose that the thermoelectric power distinguishes two competing scenarios for quantum 
phase transitions in heavy fermions : the spin-density-wave (SDW) theory and breakdown of the 
Kondo effect, fn the Kondo breakdown scenario, the Seebeck coefficient turns out to collapse from 
the temperature scale E* , associated with quantum fluctuations of the Fermi surface reconflguration. 
This feature differs radically from the physics of the SDW theory, where no reconstruction of the 
Fermi surface occurs, and can be considered as the hallmark of the Kondo breakdown theory. We 
test these ideas, upon experimental results for YbRh2Si2. 

PACS numbers: 71.27.+a, 72.15.Qm, 75.20.Hr, 75.30.Mb 



A quantum transition from a light metal to a heavy 
Fermi liquid is believed to occur in a class of heavy fer- 
mion (HF) compounds as a result of the competition 
between Rudderman-Kittel-Kasuya-Yoshida interactions 
and the formation and condensation of Kondo singlets 
[l||. The Kondo breakdown theory [2r4] offers a mini- 
mal model describing this transition. It differs from the 
spin-density- wave (SDW) theory fa, 16] in respect that the 
whole heavy Fermi surface is destabilized at the quantum 
critical point (QCP) and a reconstruction of the Fermi 
surface is inevitable, sometimes identified with an orbital 
selective Mott transition. As a result, the Kondo break- 
down QCP is described by critical fluctuations with the 
dynamical exponent z = 3, associated with breakdown 
of the Kondo effect d Q, while the SDW QCP is ex- 
plained within z = 2 related with antiferroniagnetic spin 
fluctuations jsl, ISy. 

Until now, the thermopower has been studied mainly 
in the heavy Fermi liquid phase, with a special focus on 
its sign and its ratio with the Sommerfeld coefflcient. In 
heavy fermions, the sign of the thermopower is determi- 
ned by the position of the Kondo resonance with respect 
to the Fermi surface. The strong mass renormalization is 
proven not to affect the quasi-universal ratio to the Som- 
merfeld coefficient, which remains almost the same as in 
conventional metals [7|, IsJ . Little is known about the See- 
beck coefficient close to the QCP. Preliminary studies for 
CeCu6-a;Aua; [9] and Ce(Nii_a;Pda;)2Ge2 ^ show that 
the presence of a QCP modifies low temperature depen- 
dence of the Seebeck coefficient. Two recent studies un- 
der magnetic fields show some striking similarity between 
thermoelectric effects in CeCoIus [11^ and URu2Si2[l3- 
In particular, both systems exhibit a pronounced aniso- 
tropy in their thermoelectric response. Lastly, a recent ex- 
periment on YbRh2Si2 under weak magnetic fields shows 
some drastic variations of the magnitude of the Seebeck 
coefficient in both sides of the QCP llSll . Even fewer theo- 
retical studies are available [71, [iJ, ll5| . In the case of the 



SDW the authors of Ref. [16| argued that the Seebeck 
coefficient divided by temperature has the same varia- 
tion in temperature as the Sommerfeld coefficient and 
the quasi-universal ratio is preserved at the SDW QCP. 

In this study we show that the thermoelectric power 
can be regarded as one of the hallmarks for the Kondo 
breakdown quantum criticality, providing a careful fit to 
the recent experimental observations on YbRh2Si2 |13| . 
Precisely, the abrupt collapse from the temperature iden- 
tified with E* in the Kondo breakdown scenario [3|, \^ 
turns out to be the unique prediction from the Kondo 
breakdown theory beyond the SDW [5, 6] and local quan- 
tum critical [17| scenarios. 

We start from the U(l) slave-boson representation of 
the Anderson lattice model in the large-?/ limit 
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Here, Ci^ and dio- — blfia- are conduction electron with 
a chemical potential /i and localized electron with an 
energy level ey, where bi and fia- are called holon and 
spinon, representing hybridization and spin fiuctuations, 
respectively. The spin-exchange term for the localized or- 
bital is introduced as a competitor of the hybridization 
term, and decomposed via exchange hopping processes of 
spinous, where Xij is a hopping parameter for the decom- 
position. Xi is a Lagrange multiplier field to impose the 
single occupancy constraint 6|6i -I- fj^fia — N/2, where 
N is the number of fermion flavors with a — 1, ..., N. 
Performing the saddle-point approximation of bi — >■ 5, 



Xij — >■ X^ s-nd i\i -^ A, one finds an orbital selective 
Mott transition as Kondo breakdown at J w Tr-. For 
{hi) =0, a spin- liquid Mott insulator arises with a small 
area of the Fermi surface in J > Tx while for {hi) 7^ a 
heavy Fermi liquid obtains with a large Fermi surface in 

Tk > J [i%. Here, Tk = Dexp('j^^^ is the single- 
ion Kondo temperature, where pc ~ {2D)~^ is the density 
of states for conduction electrons with the half bandwidth 
D. Reconstruction of the Fermi surface occurs at J ~ Tk- 

The fluctuation-corrections are taken into account in 
the Eliashberg framework [3, 4j . The main physics is that 
the Kondo breakdown QCP is multi-scale. The dynamics 
of the hybridization fluctuations is described by z = 3 
critical theory due to Landau damping of electron-spinon 
polarization above an intrinsic energy scale E* , while by 
z = 2 dilute Bose gas model below E* . The energy scale 
E* originates from the mismatch of the Fermi surfaces of 
the conduction electrons and spinous, one of the central 
aspects in the Kondo breakdown scenario. Physically, one 
may understand that quantum fluctuations of the Fermi 
surface reconfiguration start to be frozen at T ss E* , 
thus the conduction electron's Fermi surface dynamically 
decouples from the spinon's one below E* . We show that 
the Seebeck coefficient collapses at E* , associated with 
the Fermi surface reconstruction. 

The thermoelectric power can be deduced from the fol- 
lowing transport equations 
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defining three transport coefficients of a, p, and k, where 
three species of conduction electrons, holons, and spi- 
nons are taken into account for each response function, 
denoted by c, 6, and /, respectively. Jei(th) is the elec- 
tric (thermal) current, and E, e, /i, T are an external 
electric field, internal one, chemical potential and tempe- 
rature, respectively. The internal electric field e emulates 
the "back- flow" condition J^i -h J^; = fl3|, central to 
transport in gauge theories. Taking /ic = M/ ~ A^fc with 
an open-circuit boundary condition J^j — J^i — 0, we flnd 
the total thermopower conductivity pt{T) and electrical 
onetTt(r) 



Pt{T)^p,{T) + 
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Then, the physical Seebeck coefficient is given by St{T) 
-pt{T)/at{T). 



Around the QCP, the Seebeck coefficient can be sim- 
plified as follows 

c irr\ ~ Pc{T)+Pf{T) p^{T) 

^HF[l)'=s-—77fr——77p:, ^qc,sl(i) ga - , . ,(pJ 
ac{l) + crf{l) (Jc(l) 

based on the fact that Uf{T) <^ (^b{T) -> 00 in the HF 
phase and cr/(T) ;» <^b{T) — > in both the quantum 
critical regime and spin liquid state. 

Evaluating the thermopower conductivity in the one- 
loop level 
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where Gc(f) (fc, 1^ + i5) is the full Green's function in the 

Eliashberg approximation and Vp is the Fermi velocity 
of conduction electrons (spinous), we find 
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are wave- function and dispersion renormalization [19|, 

respectively, and '&Yjc,(f){kp ,ijj) is the real part of the 
electron (spinon) self-energy. This expression is basically 
the same as the standard representation [7|,llJ,|l5|, where 
the derivative of the density of states with respect to fre- 
quency is replaced with that of the dispersion renorma- 
lization function with respect to momentum. Thus, the 
sign of the Seebeck coefficient is given to be negative 
when the Kondo resonance lies below the Fermi energy, 
exactly the Yb case. 

An important quantity is the quasiparticle weight 
Zc(/)(w), resulting from the linear dependence of fre- 
quency in the thermopower expression. A singular loga- 
rithmic temperature dependence is revealed in the quan- 
tum critical regime, typical of the z = 3 quantum criti- 
cality in three dimensions beyond the SDW theory [16] . 

The self-energy correction can be found within the 
Eliashberg theory, where quantum corrections are self- 
consistently introduced in the one-loop level but vertex 
corrections are neglected [3|, |j|. The dispersion renor- 
malization function is obtained from the renormalized 
band. In the HF phase the spinon and conduction elec- 
tron bands are hybridized. Resorting to these two renor- 
malization functions, we can provide a fit for the recent 
experimental data for YbRh2Si2 [13]. The main feature 



of the experimental results is that the Seebeck coefficient 
divided by temperature is found to increase logarithmi- 
cally in temperature above a certain energy scale Tj^i 
and drops down abruptly in both the quantum critical 
regime and antiferromagnetic phase while it saturates to 
a constant in the heavy Fermi liquid. The Kondo break- 
down QCP scenario reproduces this main feature almost 
exactly as shown in Fig. 1. The logT divergence is explai- 
ned by the z — i quantum critical regime of the theory. 
In particular, the collapse at Tm is attributed to suppres- 
sion of quantum fluctuations of Fermi surfaces at E* , thus 
Tm = E* in the Kondo breakdown theory, where loca- 
lized f-electrons cannot participate in carrying entropy 
due to the Fermi surface decoupling. In this regime the 
question whether S/T changes sign at very low tempera- 
tures is still experimentally under scrutiny. In the Kondo 
breakdown theory, the negative sign of the Seebeck co- 
efficient is attributed to the presence of the fluctuating 
hybridization between the f-spinons and the conduction 
electrons, yielding to the formation of the Kondo reso- 
nance below the Fermi level. Below E* however, the spi- 
non Fermi surface decouples from the conduction electron 
one, and S/T saturates to the value determined solely by 
the light conduction electrons. In the heavy Fermi liquid 
the Kondo breakdown theory shows the saturation asso- 
ciated with Fermi liquid physics, but without changing 
sign. 

To understand better the quantum critical region, we 
approximate the self-energy and obtain the typical z = 3 
form analytically in three dimensions. 
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where D is the half-bandwidth and nii, — {2NV'^ pc)~^ is 
the band mass of holons. Then, the two renormalization 
functions are given by 
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displaying the log T-dependence and negative sign due to 
the f-resonance below the Fermi energy. 

Lastly, we show that in the HF regime {V > Vc), the 
saturation value of the Seebeck coefficient divided by the 

temperature a = "^^ ' is proportional to B — B*, 

where B* is slightly displaced from the critical field Be 
of the antiferromagnetic QCP in YbRh2Si2. Resorting 
to the condensation amplitude b^ « C l""^., ( -^ J (Jk — 

J^), where C is a positive numerical constant from the 
mean-field analysis, we obtain 
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Figure 1: (Color online) Fitting to the experimental data 
of Ref. rial based on the analytic expression of Eq. (0. The 
values of the parameters are the half bandwidth of the conduc- 
tion electrons D — 990K [20], of the spinons Df = IK, and 
E* = 40mA'. The upper left quadrant represents the experi- 
mental data. The upper right quadrant fits the QC regime. 
An overall adjusting coefficient A = 0.33 has been used. The 
holon mass is m = 6.0 x 10"*-?^. The lower left quadrant cor- 
responds to a field of lOOmT. A coefiicient A = 0.3 has been 
used with holon mass m — 5.5 x 10~*_R'. The lower right qua- 
drant corresponds to a field of IT. A coefficient A = 0.0135 
has been used with holon mass m — —4.1 x lO"*-/^. The logT 
behavior of the z — 3 QC regime is correctly reproduced while 
the scale E* corresponding to the reconfiguration of the Fermi 
surface is found to be a bit small. 



where J§ 
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J^ . Assuming 



Jk oc B valid near the QCP, we conclude a oc B 
consistent with the experimental data of Ref. [13]. 

The weak point of the Kondo breakdown theory is 
the treatment of the anti-ferromagnetism. However, the 
thermoelectric power is insensitive against onset of anti- 
ferromagnetism as far as E* is larger than T^, Neel tem- 
perature. Since localized f-electrons cannot carry entropy 
when they are decoupled from conduction electrons in 
the quantum level, the collapse should occur from E* 
above Tn- Actually, this is observed in YbRh2Si2, where 
E* « lOOmK and T/v ~ 70mK as the maximum Neel 
temperature [13|. This collapse behavior is believed to 
survive beyond our approximation. 

This feature differs from the SDW scenario completely, 
well discussed in Ref. [2ll|. Because quantum fluctuations 
associated with the Fermi surface reconfiguration do not 
exist in the SDW theory, the Seebeck coefficient should 
saturate to a constant value below both the Neel tempe- 
rature and HF coherence one. It should display a sym- 
metric configuration for B < Be and B > Be- 

Another important result of the Kondo breakdown 
theory applied to YbRh2Si2 is that the sign of the See- 
beck coefficient becomes positive in the low temperature 
limit, changing from that in the HF phase jl3|. Since 
the spinon band is decoupled from the conduction band. 



the Kondo resonance disappears even in the dynamical 
level and contributions from only conduction electrons 
give rise to a small positive constant, reflecting the phy- 
sics of the normal Fermi liquid. This feature is difficult to 
interpret within the SDW framework without the Fermi 
surface reconstruction. 

One may ask why the Seebeck coefficient is sensitive 
to E* while other quantities such as the specific heat co- 
efficient 22| and thermal conductivity [23| do not show 
such serious dependence. The main difference between 
the thermoelectric power and others is that the thermo- 
dynamic and thermal transport coefficients can measure 
contributions of neutral spinon excitations while the See- 
beck coefficient does not. This is the reason why only the 
thermopower is sensitive on the Fermi surface reconstruc- 
tion, argued to be a fingerprint for the Kondo breakdown 
QCP. 

Another important measurement is the Hall coefficient 
which also reveals an interesting energy scale, referred as 
T* in the heavy fermion side, below which the Hall coef- 
ficient Rh shows an abrupt decrease. This abrupt change 
has been attributed to the proximity to a Kondo break- 
down QCP |24l . l25fl , whereas other approaches have stres- 
sed the complexity of this measurement and its sensitivity 
to minus fluctuations of the f-electron chemical potential 
[2d |. Since the Hall coefficient measures the Fermi sur- 
face curvature, thus sensitive to the static formation of 
the heavy fermion band, the Hall number shows its cha- 
racteristic only in the heavy fermion side while it does 
not in the localized side. On the other hand, the ther- 
moelectric power measures fluctuations of Fermi surfaces, 
showing an interesting signature in the antiferromagnetic 
side. We refer interested readers to our future publication 
for a thorough study of the Hall effect within the Kondo 
breakdown theory. 

In this paper, we argue that the thermoelectric power 
is an important measurement for quantum fluctuations 
of the Fermi surface reconfiguration, which enables us to 
discriminate the Kondo breakdown theory from the SDW 
scenario. It was demonstrated to collapse at the energy 
scale associated with Fermi surface fluctuations, identi- 
fied with E* in the Kondo breakdown theory. We show 
here, that a careful fitting of the data can be obtained 
within the Kondo breakdown theory for an energy scale 
of the order of 40mK, a bit small compared to the expe- 
rimental value. The z = 3 quantum criticality of Fermi 
surface fluctuations gives rise to the singular log T beha- 
vior in the quantum critical regime, consistent with the 
experimental data on YbRh2 Si2 . 

The local quantum critical scenario [17| does not have 
the characteristic energy scale E* of the Kondo break- 
down theory. There, the mechanism for the abrupt col- 
lapse of the Seebeck coefficient is not yet as clear as wi- 
thin the Kondo breakdown theory. It would be very inter- 
esting to see how the two theories of Kondo breakdown 
compare with respect to the fitting of the experimental 



data. 
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